Theoretical derivation of 1// noise in quantum chaos 
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It was recently conjectured that 1/f noise is a fundamental characteristic of spectral fluctuations 
in chaotic quantum systems. This conjecture is based on the behavior of the power spectrum of 
the excitation energy fluctuations, which is different for chaotic and integrable systems. Using 
random matrix theory we derive theoretical expressions that explain the power spectrum behavior 
at all frequencies. These expressions reproduce to a good approximation the power laws of type 
1/f (1//^) characteristics of chaotic (integrable) systems, observed in almost the whole frequency 
domain. Although we use random matrix theory to derive these results, they are also valid for 
semiclassical systems. 
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The study of energy level fluctuations is a basic tool in 
understanding quantum chaos. The pioneering work of 
Berry and Tabor showed that the spectral fluctuations of 
a quantum system whose classical analogue is fully inte- 
grable are well described by Poisson statistics Q . On the 
other hand, Bohigas et al. conjectured that the fluctua- 
tion properties of quantum systems which in the classical 
limit are fully chaotic coincide with those of Random Ma- 
trix Theory (RMT) 0. Recently, a different approach to 
quantum chaos has been proposed Considering the 
sequence of energy levels as a discrete time series, it has 
been conjectured that chaotic quantum systems are char- 
acterized by 1// noise, whereas integrable quantum sys- 
tems exhibit 1//^ noise. This conjecture is supported by 
numerical experiments which involve classical Random 
Matrix Ensembles (RME) and atomic nuclei 

In this Letter we present a theoretical derivation of 
1// and 1//^ noises in chaotic and regular systems, re- 
spectively. We use RMT to derive these results, but they 
are also valid for semiclassical systems. We present the 
main steps of the derivation and compare the theoretical 
results with numerical calculation for RME, an atomic 
nucleus and a quantum billiard, finding excellent agree- 
ment. 

For any quantum system the accumulated level density 
N{E) can be separated into a smooth part N{E) and a 
fluctuating part N{E). To remove the main trend defined 
by the former, the energy levels Ei are mapped onto new 
dimensionless levels = N{Ei). This transformation, 
called unfolding, allows to compare the statistical prop- 
erties of different systems or different parts of the same 
spectrum. 

The analogy between the energy spectrum and a dis- 
crete time series is established in terms of the 5q statistic. 



defined as 

<? 

i=l 

where is the i-th unfolded level and Si = ei+i — Cj. 
Once the unfolding is performed, the average nearest 
level spacing is (s) ~ 1. Note that 5q represents the 
deviation of the excitation energy of the (g -I- l)-th un- 
folded level from its mean value. Moreover, it is closely 
related to the level density fluctuations. Indeed, we can 
write 

6q = -iV(i?,+i), (2) 

if we appropriately shift the ground state energy; thus, it 
represents the accumulated level density fluctuations at 
E = Eq+i. Its power spectrum, defined as the square of 
the modulus of the Fourier transform, shows neat power 
laws both for fully chaotic and integrable systems, i.e., 



but level correlations change the exponent from a = 2 
for uncorrelated spectra to a = 1 for chaotic quantum 
systems 0. 

Notation.- We consider an interval of length L S> 1 
containing N c^i L unfolded energy levels. The fluctuat- 
ing parts of the level and accumulated level densities are 
denoted p(e) and n(e), respectively. In addition to the 
5q statistic, wc also introduce another discrete function 
Uq = n{q), obtained by sampling the continuous function 
for integer values of the energy. The Fourier transforms 
and power spectra of these functions are dcflned in the 
usual way The differences between continuous and 
discrete functions, and specially between their Fourier 
transforms, play a relevant role in the following. The no- 
tation we use for all these functions, Fourier transforms 
and power spectra is summarized in Table 
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TABLE I: Summary of the functions, Fourier transforms and 
power spectra used in this letter 



Domain 


R 


Z 


Function 


n{e) 


riq 






Fourier transform 


n(r) 


rik 






Power spectrum 


P"(t) 


Pk 




Pi 



Spectral and ensemble averages will be denoted by (•), 
or (•)^ to distinguish ensemble averages performed in dif- 
ferent RME. Here, /3 stands for the repulsion parame- 
ter characterizing the ensemble. In this work wc shall 
consider two RME: the Gaussian Orthogonal Ensemble 
(GOE) and the Gaussian Unitary Ensemble (GUE) 
Wc have /3 = 1 for GOE and /3 = 2 for GUE. 

Spectral fluctuations. - The main object of this Letter 
is to obtain explicit expressions of the average value of 
for chaotic and integrable systems. Except for integrable 
systems, one of the main features of quantum spectra is 
that successive level spacings are not independent ran- 
dom variables, but correlated quantities. This property 
makes exceedingly difficult to work directly with the dis- 
crete Sq sequence. To circumvent this difficulty we profit 
from the relationship (|2J). The statistical properties of 
n(e) are usually measured in terms of the spectral form 
factor, defined as 



Kir) 



dep{e) 



(4) 



that is, as the power spectrum of the fluctuating part 
of the energy level density. Instead of K{t), wc can use 
the power spectrum P"{t) of n(e) to analyze spectral 
fluctuations. It can be shown that for very large L values 



i 



(P"(r)) 



Kir) 



(5) 



Then, we shall (i) sample n(e) for integer values of the 
energy and study the fluctuations of the new riq function; 
(ii) map the independent variable of (an energy) to the 
independent variable of Sq (a dimensionless index). 

Quantum Chaos.- One of the most important features 
of fully chaotic system is the universal behavior of their 
spectral fluctuations. Therefore, we can use RMT to 
obtain the main properties of {P"{t)) in this kind of 
systems. Exact analytical expressions are known for 



K{t) in all RME, leading to 



, r<l, 
r, r>l, 



(6) 



47r2r2 ' 
2 - log(l + 2r) 



2 - r log 



2t + 1 
2t- 1 



r< 1, 



, r > 1. 



(7) 



For small r values, © and Q become (see Ref. 0) 



1 



T< 1, 



'1^ 2/3^V 
and for t > 1 we can approximate (^'"(t))^ by 



(P"(r)) 



1 



T> 1, 



(8) 



(9) 



which is exactly equal for GUE and a good approximation 
for GOE. 

On the other hand, using periodic orbit theory and 
semiclassical mechanics it is possible to calculate K{t) for 
chaotic systems. The semiclassical expression essentially 
coincides with the results of RMT for Tmin ^ t <^ th, 
where Tmin is the period of the shortest periodic orbit, 
and th — h/ (s) is the Heisenberg time of the system, 
related to the time a wave packet takes to explore the 
complete phase space of the system 0, . The expres- 
sions derived in the following are then directly applicable 
to generic chaotic quantum systems for t between this 
two values. 

The next step is to relate the spectral fluctuations of 
riq to those of the continuous function n(e), as given by 
their power spectra. With the usual definitions for the 
continuous and discrete Fourier transforms we have 



1 



V 5 
N ^ 

q— — OQ 



k 



(10) 



Therefore, to relate and P"{t) wc need the precise 
knowledge of n(r) for all r. Since RMT only provides 
the mean value of its squared modulus, we introduce 
two simplifying assumptions: (a) for times r — k/N, 
mk/N)\ = ^/WJP^Hk/N)); (b) the phases of n{k/N) 
are random variables uniformly distributed in the inter- 
val [0, 27r). These assumptions, specially (a), are reason- 
able as far as we are only interested in the average values 
of P"(t) and PJ}. With these simplifications and using 
eq. © we obtain the following result for 1, 



(P;!') = (P"(fc/A^))-f (P"(l-fc/iV))- 

1 



47r2(Ar- fc)2 



4smM- 



47r2fc2 
fc = l,2,. 



,7V- 1. 



(11) 
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Finally, from the relationship between the variances of 
6q and n(e) for chaotic systems |loj we obtain an analyt- 
ical expression for (P^) in terms of {P^{k/N)); the final 
result is 

(P^) = {PI}) - 1^, fc=l,2,... ,iV-l, iV»l. (12) 

Eq. ifT^ . together with and (7J) or gives ex- 
plicit expressions of (P^) for GOE and GUE. For generic 
chaotic systems we can apply GOE or GUE results de- 
pending on their symmetries. 

Integrable systems.- A similar calculation can be per- 
formed for integrable systems. In this case K{t) = 1 [llj, 

and instead of (El or Q we have P"(t) = — ^ m for 

all T. Using obtained with the same assumptions 

(a) and (b), and the fact that the variances of Sq and n{e) 
are equal for integrable systems, we get 



(no 



1 



4smM- 



> 1. 



(13) 



In order to check the assumptions (a) and (b) intro- 
duced above, we may now give an alternative proof of 
eq. H13|l without using these assumptions. For inte- 
grable systems, the spacings sequence Si can be con- 
sidered as a sequence of independent random variables 
with Poisson distribution Q. Therefore, the 6q statis- 
tic is a sum of such random variables Sq = '^^^I'lVi, 
with Wi = Si— < s >. Additionally, in this case we 

|2\ 



have (wk) = and {^\wk\ / ~ where Wk is the Fourier 
transform of Wj. Then, we can write 



1 



Af-l 



q 



/ ,Wk } 

k=0 m=l 



2'nikm 
N 



(14) 



This series can be added up considering that {wk) — 0; 
we obtain for A^ 1 



-ink 



1 X ^ lUUkC JV 

7n ^ 



fc=o 2 sin 



nk 

TV 



2'iTikq 
N . 



(15) 



and therefore 



5k = 



iWke 



N 



, Trfc 
2sin(- 



A^> 1. 



(16) 



Consequently, taking into account that for a Poissonian 
sequence (|wfc|^ ) = 1, wc can write 



{Pi) 



Sk 



1 



4sm^ ( ^ 



Af > 1. 



(17) 



The coincidence of eqs. (|13|l and H17() shows that the 
assumptions (a) and (b) introduced to derive H12() and 
are sound approximations. 
1// and 1/p noises.- Note that eqs. (O and ifT^ 
are generic results, valid for every chaotic and integrable 
quantum system, respectively. Nevertheless, there can 
be some deviations at the lower frequencies due to the 
behavior of the shortest periodic orbits. When k <^ N 
the first term of eq. (|ll|l becomes dominant and using 
eq. we can write, for chaotic systems. 



(Pi) - 



N 



2f3'!T^k 



, fc < A, TV> 1. 



Similarly, for integrable systems, eq. H13|l becomes 



(Pi) = 



(18) 



(19) 



These expressions show that for small frequencies, the 
excitation energy fluctuations exhibit 1// noise in chaotic 
systems and 1//^ noise in integrable systems. As we 
shall see below, these power laws are also approximately 
valid through almost the whole frequency domain, due to 
partial cancellation of higher order terms in eqs. (|12|l and 
Only near k = N/2, the so called Nyquist frequency 
the effect of these terms becomes appreciable. 




log(k) 

FIG. 1: Theoretical power spectrum of the Sq function for 
GUE and integrable systems (solid lines), compared to the 
results calculated numerically using 30 GUE matrices of di- 
mension N = 10^ (circles) and 30 Poisson level sequences of 
length N = 10^ (triangles). 



To test all these theoretical expressions we have com- 
pared their predictions to numerical results obtained for 
different ensembles and systems. Fig. ^ displays the 
theoretical values of (Pk) for GUE and integrable sys- 
tems, as given by (fT^ and (fT!^ . together with the nu- 
merical results for GUE matrices and Poisson level se- 
quences published in In order to enlarge the high 
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frequency region, where the numerical results show small 
deviations from the 1//" power law behavior, an upper 
right panel is added to the figure. The agreement be- 
tween the theoretical and numerical results is excellent 
at all frequencies (note that there are no free parameters 
in the analytical result). The theoretical curve describes 
perfectly the power laws, characteristic of small and in- 
termediate frequencies, as well as the deviations observed 
at the highest frequencies. 
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FIG. 2: Numerical average power spectrum of the 5q function 
for ''^Na, calculated using 25 sets of 256 consecutive levels 
from the high level density region, compared to the parameter 
free theoretical values (solid line) for GOE. 
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FIG. 3: Numerical average power spectrum of Sq for a rectan- 
gular billiard, calculated using 25 sets of 256 consecutive lev- 
els, compared to the parameter free theoretical values (solid 
line) for integrable systems. 

We have also compared our predictions to the power 
spectra of Sq for two physical systems: an atomic nucleus 
(chaotic) and a rectangular billiard (integrable). 



In the first case we have performed a shell model cal- 
culation for '^*Na using an adequate realistic interaction. 
The Hamiltonian matrices for different angular momenta, 
parity and isospin were fully diagonalized. Then, 25 sets 
of 256 consecutive high energy levels of the same J^T 
were selected, and the average power spectrum of the 
6q function was calculated. Fig. |21 shows the result of 
this calculation together with the theoretical values of 
eq. H12|) . An excellent agreement between the theoreti- 
cal and numerical results is obtained through the whole 
frequency interval. 

As an example of integrable system we have chosen 
a rectangular billiard with sides of length a = v'A and 
b = l/yX, with A = (75 + l)/2; this geometry gives 
rise to an irrational ratio a/b = A, and thus there are 
no degeneracies in the spectrum. We have calculated the 
spectrum and selected 25 sets of 256 consecutive very 
high energy levels in order to avoid, as far as possible, 
the influence of short periodic orbits. Fig. O shows the 
results for the average power spectrum of Sq and the the- 
oretical values given by eq. ((T^ . As in previous cases, we 
can see that the agreement between the theoretical and 
numerical results is very good in the whole frequency do- 
main. 

In summary, we have derived theoretical expressions 
for the power spectrum of the Sq function both for regu- 
lar and chaotic quantum systems. These expressions are 
universal and do not contain any free parameter. We 
have compared our theoretical predictions with numeri- 
cal results for RME, a rectangular billiard and an atomic 
nucleus, obtaining an excellent agreement for all these 
systems. The theory reproduces the power laws of type 
1// for chaotic systems, and 1//^ for regular ones, ob- 
served in the power spectrum of the excitation energy 
fluctuations up to frequencies very close to the Nyquist 
limit. 
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